View(final)
final$labels = paste(final$quant,   "alpha = ", round(final$alpha, 2))
View(final)
final$labels = paste(final$quant,   "(alpha = ", round(final$alpha, 2), ")")
View(final)
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
final$labels = paste(final$quant,   expression(alpha), round(final$alpha, 2), ")")
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
bquote(.(final$quant) ~~~ alpha)
test = bquote(.(final$quant) ~~~ alpha)
test = list(bquote(.(final$quant) ~~~ alpha))
test = list(bquote(alpha == .(final$alpha)))
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
View(final)
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
View(labs)
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
test = list(bquote(alpha == .(final$alpha)))
test
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
test
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
unique(final$quant)
test
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
alpha = unique(final$alpha)
alpha
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
library(here)
dir = here()
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
gA = ggplotGrob(firm_plot)
gB =  ggplotGrob(quant_plot)
setwd(dir)
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_v_firm_size.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/gem_power_law.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
View(merge)
View(final)
wd = gsub("Figures", "Empirical Data/GEM", dir)
library(data.table)
library(ggplot2)
library(here)
library(gridExtra)
library(grid)
library(gtable)
library(magrittr)
dir = here()
# Firm size
wd = gsub("Figures", "Empirical Data/GEM", dir)
setwd(wd)
GEM = fread("energy_v_firm_size.csv")
GEM$year_country = paste(GEM$code, GEM$year)
wd
GEM = fread("energy_v_firm.csv")
x = log(GEM$energy)
y = log(GEM$firm_mean)
r = lm(y ~ x)
r.squared = summary(r)$r.squared
r.squared = paste("R^2 == ", round(r.squared,2) )
x =  exp(seq(log(1), log(10^5), length.out = 40 ))
x = log(x)
x = data.frame(x)
p = predict(r, x,  interval="confidence", level = 0.99)
p = exp(p)
p = data.frame(x = exp(x), p)
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_v_firm_size.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_v_firm_size.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
View(GEM)
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
View(GEM)
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
View(GEM)
quant.func = function(x){
c(min(x), max(x))
}
quant.func = function(x){
c(min(x), max(x)) %>% mean()
}
x = 1:10
quant.func(x)
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
GEM[, mean(energy), by = energy_quintile]
GEM[, mean(energy), by = energy_quantile]
GEM[, as.list(min(energy),max(energy)), by = energy_quantile]
test = GEM[, as.list(min(energy), max(energy)), by = energy_quantile]
View(test)
test = GEM[, max(energy), by = energy_quantile]
test2 = GEM[, min(energy), by = energy_quantile]
look = data.table(test, test2)
View(look)
test = GEM[, mean(energy), by = energy_quantile]
plot(test$energy_quantile, test$V1)
test = GEM[, mean(energy), by = energy_quantile]
plot(test$energy_quantile, test$V1)
View(test)
plot(test$V1)
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
quantile_mean = GEM[, .(mean = mean(energy)), by = energy_quantile]
View(quantile_mean)
bar_plot = ggplot() +
geom_bar(data = quantile_mean, (x = energy_quantile, y = mean))
bar_plot = ggplot() +
geom_bar(data = quantile_mean, aes(x = energy_quantile, y = mean) )
bar_plot
bar_plot = ggplot() +
geom_bar(data = quantile_mean, aes(x = energy_quantile, y = mean), stat="identity" )
bar_plot
bar_plot = ggplot() +
geom_bar(data = quantile_mean, aes(x = energy_quantile, y = mean, col = energy_quantile), stat="identity" )
bar_plot
bar_plot = ggplot() +
geom_bar(data = quantile_mean, aes(x = energy_quantile, y = mean, fill = energy_quantile), stat="identity" )
bar_plot
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
inset = ggplotGrob(bar_plot)
piketty_hist = quant_plot + annotation_custom(grob = inset, xmin = 1.2, xmax = 10, ymin = -5, ymax = -3)
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/energy_firm_power.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/case_pay_shape.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/case_pay_shape.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/case_shape.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/case_pay_shape.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/case_shape.R')
library(ggplot2)
library(gridExtra)
library(data.table)
library(hmod)
library(magrittr)
library(scales)
library(here)
text.size = 10
dir = here()
wd = gsub("Figures", "Hierarchy Model/results", dir)
setwd(wd)
mod = fread("power_concentration.csv")
names(mod) = c("firm_size", "power_gini")
wd = gsub("Figures", "Empirical Data/GEM", dir)
setwd(wd)
setwd(wd)
energy = fread("energy_v_firm.csv")
source('~/Desktop/origin_inequality/Supplementary Material/Figures/power_concentration.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/power_concentration.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/lenski_plot.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/lenski_plot.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/lenski_plot.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/lenski_plot.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/lenski_plot.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/lenski_plot.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/lenski_plot.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/power_concentration.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
n = 10000
alpha = 2.1
large_firm_sample = rpld(10^7, 1, alpha, 2.3*10^4)
high_mean_firm = mean(large_firm_sample)
high_energy = e_predict(high_mean_firm) %>% round()
high_energy_label = paste(  "Industrial Society \nEnergy Use per Capita ~ ", high_energy, "GJ per person")
firm_sim = sample(large_firm_sample, n)
base_sim = base_fit(a, b, firm_sim)
r_sim = r_sim_2(fit$V1, fit$r, firm_sim)
base_pay_sim = bp_sim(fit$base_pay, n)
mod_high = model(a, b, base_sim, firm_sim, base_pay_sim, r_sim, sigma, power = T) %>% data.table()
high_power_gini = fastgini(mod_high$power)
high_power_gini = fastgini(mod_high$power) %>% round(.,2)
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
high_energy = e_predict(high_mean_firm)
high_energy = round(high_energy/100) * 100
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
firm_sim = sample(large_firm_sample, 10^6)
base_sim = base_fit(a, b, firm_sim)
r_sim = r_sim_2(fit$V1, fit$r, firm_sim)
base_pay_sim = bp_sim(fit$base_pay,10^6)
mod_high = model(a, b, base_sim, firm_sim, base_pay_sim, r_sim, sigma, power = T) %>% data.table()
high_power_gini = fastgini(mod_high$power) %>% round(.,2)
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_landscape_growth.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_hierarchy_shape.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_hierarchy_shape.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_hierarchy_shape.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/mod_hierarchy_shape.R')
power = c(1,1,2)
g = Gini(g)
source('~/Desktop/origin_inequality/Supplementary Material/Figures/power_gini_example.R')
g = Gini(power)
?Gini
g = Gini(power, corr = T)
power = c(1,1,3)
g = Gini(power, corr = T)
power = c(rep(1, 4), 3, 3, 7)
g = Gini(power, corr = T)
power = c(rep(1, 8), rep(3, 4) , 7, 7, 15)
g = Gini(power, corr = T)
power = c(rep(1, 16), rep(3, 8), rep(7, 4) , 15, 15, 31)
g = Gini(power, corr = T)
power = c(1,1,3)
g = Gini(power, corr = T)
power = c(1,1,3)
rep(power, 10^4)
powe = rep(power, 10^4)
g = Gini(power, corr = T)
power = c(1,1,3)
power = rep(power, 10^4)
g = Gini(power, corr = T)
power = c(1,1,3)
g = Gini(power, corr = T)
power = c(1,1,3)
g = Gini(power)
power = c(1,1,3)
power = rep(power, 10^4)
g = Gini(power)
power = c(rep(1, 4), 3, 3, 7)
g = Gini(power)
power = c(rep(1, 4), 3, 3, 7)
g = Gini(power)
power = c(rep(1, 8), rep(3, 4) , 7, 7, 15)
g = Gini(power)
power = c(rep(1, 16), rep(3, 8), rep(7, 4) , 15, 15, 31)
g = Gini(power)
power = c(1,1,3)
power = rep(power, 10^4)
g = Gini(power)
power = c(rep(1, 4), 3, 3, 7)
g = Gini(power)
power = c(rep(1, 8), rep(3, 4) , 7, 7, 15)
g = Gini(power)
power = c(rep(1, 16), rep(3, 8), rep(7, 4) , 15, 15, 31)
g = Gini(power)
power = c(rep(1, 32), rep(3, 16), rep(7, 8), rep(15, 4), 31, 31, 63)
g = Gini(power)
source('~/Desktop/origin_inequality/Supplementary Material/Figures/power_gini_example.R')
hist(power)
library(ggplot2)
library(readr)
library(gridExtra)
library(grid)
library(gtable)
library(data.table)
library(hmod)
library(snow)
library(doSNOW)
library(foreach)
library(tcltk)
library(here)
n.iteration = 10000
dir = here()
wd = gsub("Figures", "Hierarchy Model/data", dir)
setwd(wd)
s = fread("s_empirical.txt")
g = fread("g_empirical.txt")
source('~/Desktop/origin_inequality/Supplementary Material/Figures/case_parameters.R')
library(ggplot2)
library(boot)
library(magrittr)
library(data.table)
library(dplyr)
library(grid)
library(ineq)
library(here)
text.size = 10
dir = here()
wd = gsub("Figures", "Empirical Data/Case Studies/Number of Subordinates", dir)
setwd(wd)
result = fread("Subordinates Results.csv")
# plot
x = log(result$N.sub)
y = log(result$P.norm)
r = lm(y ~ x)
x = seq(log(0.1), log(10^4.5), length.out = 100) %>%  data.frame  %>%  `colnames<-`(c("x"))
exp(-0.09)
